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Zitterbewegung and its effects on electrons in semiconductors 
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An analogy between the band structure of narrow gap semiconductors and the Dirac equation 
for relativistic electrons in vacuum is used to demonstrate that semiconductor electrons experience 
a Zitterbewegung (trembling motion). Its frequency is ujz ~ £ g /h and its amplitude is Xz, where 
Xz = h/rriQU corresponds to the Compton wavelength in vacuum (E g is the energy gap, mo is the 
effective mass and u ta 1.3 x 10 8 cm/sec). Once the electrons are described by a two-component 
spinor for a specific energy band there is no Zitterbewegung but the electrons should be treated as 
extended objects of size Xz- Possible consequences of the above predictions are indicated. 
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It was noted in the past that the E(k) relation be- 
tween the energy E and the wavenumber k for electrons 
in narrow-gap semiconductors (NGS) is analogous to that 
for relativistic electrons in vacuum [1-4]. The analogy 
was also shown to hold for the presence of external fields 
which was experimentally confirmed on InSb [5]. This 
" semirelativity in semiconductors" is valid for time de- 
pendent phenomena as well, so that the cyclotron res- 
onance of conduction electrons in InSb could be inter- 
preted in terms of the time dilatation between a moving 
electron and an observer [5]. The semirelativistic phe- 
nomena appear at electron velocities of 10 7 — 10 8 cm/sec, 
much lower than the light velocity c. The reason for this 
is that the maximum velocity u in semiconductors, which 
plays the role of c in vacuum, is about 10 s cm/sec. 

Until present the semirelativistic considerations for 
semiconductors were concerned with phenomena related 
mostly to classical mechanics. The purpose of this contri- 
bution is to investigate the quantum domain described by 
the Hamiltonian for energy bands in NGS, which bears 
close similarity to the Hamiltonian for relativistic elec- 
trons in vacuum. The effects we predict should be much 
more readily observable in NGS than in vacuum so this 
investigation is of interest not only for the solid state 
physics but also for the high energy physics. 

We begin with the k • p approach to InSb-type III-V 
semiconductor compounds, first written down by Kane 
[6]. Taking the limit of large spin-orbit energy A, the 
resulting dispersion relation for the conduction and the 
light-hole bands is E = ±E p , where 



E p = 



V 
'2m* 



1/2 



(1) 



Here £ g is the energy gap and m^ is the effective mass 
at the band edge. This expression is identical with the 
relativistic relation for electrons in vacuum, with the cor- 
respondence: Eg — + 2tooc 2 and — > mo- The elec- 
tron velocity v in the conduction band described by 
Eq.(l) reaches the saturation value as p increases. This 
can be seen directly by calculating Vi = dE p /dpi and 



taking the limit of large pi, or by using the analogy: 
c = (2m c 2 /2m ) 1/2 -► {£ g /2m* ) 1 / 2 = u. Taking the 
experimental parameters £ g and mj we calculate very 
similar value of u sw 1.3 x 10 8 cm/sec for different III-V 
compounds. 

Now we define an important quantity 



Az = 



m^u 



(2) 



which we call the length of Zitterbewegung for reasons 
given below. Here we note that it corresponds to the 
Compton wavelength A c = h/mac for electrons in vac- 
uum. Let us suppose that we confine an electron to the 
dimensions Az w Xz/2. Then the uncertainty of mo- 
mentum is Ap z > h/ Az and the resulting uncertainty of 
energy AE « (Ap z ) 2 /2m,Q becomes 



AE > 2m* u 2 = £ g 



(3) 



Thus the electron confined to Az w Xz/2 has the un- 
certainty of energy larger than the gap, so that it " does 
not know" whether it belongs to the conduction or to the 
valence band. 

Next we consider the band Hamiltonian for NGS. The 
momentum p must be kept in operator form, so this 
Hamiltonian is not identical with the k • p scheme of 
Kane's. Still, it is derived within the same model in- 
cluding Tq (conduction), T$ (light and heavy hole) and 
T-? (split-off) bands and it represents an 8x8 operator 
matrix [7]. We assume, as before, A>? 9 and omit the 
free electron terms since they are negligible for NGS. The 
resulting 6x6 Hamiltonian has ±£ g /2 terms on the diag- 
onal and linear pi terms off the diagonal, just like in the 
Dirac equation for free electrons. However, the three 6x6 
matrices on multiplying the momentum components Pi 
do not have the properties of 4x4 Dirac matrices, which 
considerably complicates calculations. For this reason, 
with only a slight loss of generality, we take p z and 
Px — Py — 0. In 0:3 matrix, two rows and columns corre- 
sponding to the heavy holes contain only zeros and they 
can be omitted. The remaining Hamiltonian for the con- 
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duction and the light hole bands reads 



H = ua 3 p z + -£ g [3 , 



(4) 



where 0:3 and (3 are the well known 4x4 Dirac matrices [8] . 
The Hamiltonian (4) has the form appearing in the Dirac 
equation and in the following we can use the procedures 
of relativistic quantum mechanics (RQM) . 

The electron velocity is z — (l/ih)[z, H] = ua^. The 
eigenvalues of are ±1, so that the eigenvalues of i 
are, paradoxically, ±u. In order to determine 0:3 (t) we 
calculate d 3 (t) by commuting d 3 with H and integrating 
the result with respect to time. This gives z(t) and we 
calculate z(t) integrating again. The final result is (cf. 
Refs [9-11]) 

, . . . u 2 p z . ihu ,—2iE v t. 
z(t)=z(0) + ^t + A o —e X p{-—2-) , (5) 

where Aq = 013(0) — up z /E p , and E p is given by Eq.(l) 
with p — p z . The first two terms of Eq.(5) represent 
the classical electron motion, u 2 p z /E p being the clas- 
sical constant velocity. The third term describes time 
dependent oscillations with a frequency of luz ~ £g/fi. 
Since A w 1, the amplitude of oscillations is hu/2E p w 
h/2m Q 'u = Xz/2. In RQM the analogous oscillations are 
called Zitterbewegung (ZB, trembling motion), which ex- 
plains the name given above to Xz ■ We note that the ZB 
does not obey Newton's first law since we have a non- 
constant velocity without a force. In RQM it is demon- 
strated that the ZB is a result of interference between 
states of positive and negative electron energies. 

In order to investigate the case of a definite sign of 
energy, we expand the multiband Hamiltonian of Ref [7] 
into an effective one-band Hamiltonian for the Tq con- 
duction band (here we assume quite generally p ^ 0). 
We do this by solving the initial set of six equations for 
the envelope functions (for A>^ s ) by substitution, i.e. 
expressing the valence functions fs-.-fe by the conduction 
functions fa and fa for spin-up and spin-down states, re- 
spectively. When doing this we account for the changed 
normalization condition (see [12, 13]). The final result is 



H = -S-+V+ 

2 2m* 



+ 



9 A*s 



1 , P 



9 



E g 2m 



4£ g m 



-a • ( W x P) 



A8 g m* 



V 2 V 



(6) 



where P = p + eA is the canonical momentum, V is an 
external potential, <x is the Pauli spin operator, /x_b is the 
Bohr magneton, B is an external magnetic field, and g* 
is the Lande spin factor. The first four terms represent 
the effective Pauli equation, the fifth is a nonparabolic 
correction to the energy, the sixth is an effective spin- 
orbit interaction, and the last one is the effective Darwin 
term [14]. The Hamiltonian (6) is in relativistic analogy 



to v 2 /c 2 expansion of the Dirac equation, known from 
RQM (see [15]). In our case this expansion is of the 
order v 2 /u 2 . 

The effective Darwin term can be interpreted in terms 
of ZB. If an electron oscillates around the position r, 
the potential energy can be expanded as follows: V(r + 
Ar) w V(r) + Ar • W + (Ar • V)(Ar • V)V/2. On 
average the odd powers of Ar vanish and the last term is 
(Ar) 2 V 2 V r /6. It becomes equal to the effective Darwin 
term in Eq.(6) if (Ar) 2 = (3/4)A 2 ? , which corresponds 
to Az — Xz/2. This is in agreement with our previous 
arguments. However, once the electron is described by 
a two-component spinor for a specific energy band, as 
in Eq.(6), the quantity Ar is not the amplitude of ZB 
anymore. We discuss this below. 

The four-component wave functions resulting from 
Eq.(4) can be transformed exactly into two-component 
functions for positive (or negative) electron energies. 
For free electrons this is done by applying the Foldy- 
Wouthuysen (FW) unitary transformation [16]. In our 
one-dimensional treatment the transformation reads 



£lS ^ E p + f3H 



[2E p (E p + m*u 2 )] 1/2 



(7) 



The transformed wave functions W(z) are obtained from 
the initial functions ^>(z') as follows 



*'(*) = J K{z,z')V{z')dz' , 



(8) 



where 



2E p i 



1/2 



1 . Jff 
1 + I x 



+ m*u 2 J 2 1' E p , 



exp 



ip' 3 (z-z') 



h 



d P ' 3 



(9) 



The kernel K(z, z') is not a point transformation. Sup- 
pose we are interested in the eigenfunction of position 
z. In the initial representation this function is ^f(z') — 
6(z' — zq), and it follows from Eq.(8) that the trans- 
formed function is K(z, z ). The integral of K(z, zq) over 
Z = z — z axis is unity. To get an idea of the extension of 
K(z, zq), we calculate its second moment and after some 
manipulations we obtain (see Ref. [17]) 



/ 



Z 2 K(z,z a )dZ = -A| 



(10) 



Thus the extension of the transformed eigenfunction of 
position is \z — z \ ~ Xz/2. One can show that in the 
transformed state there is no Zitterbewegung since the 
FW transformation eliminates the negative energy com- 
ponents of the wave functions [11]. All in all, following 
the interpretation established in RQM, we are confronted 
with the following choice. 1) We deal with a point-like 
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electron described by a four-component function, which 
experiences the ZB with the radius of Xz/2. 2) The elec- 
tron is described by a two-component spinor for either 
positive or negative energy and there is no ZB, but the 
electron is smeared to an object of the radius Az/2. It is 
this smearing that enters into the effective Darwin term 
in Eq.(6), as this equation describes electrons using two- 
component spinors for positive energies. In other words, 
the Darwin term is a direct consequence of smearing. 

The above result has far reaching consequences. In 
most cases electrons in a specific band are described 
by effective band parameters obtained with the use of 
the Luttingcr-Kohn (LK) transformation [18]. In fact, 
our Eq.(6) can also be obtained by this method. The 
LK transformation separates energy bands and it corre- 
sponds directly to the Foldy-Wouthuyscn transformation 
in RQM. For the FW transformation the expansion pa- 
rameter is Sfw = p/2mac = k\ c /2, while for the LK 
transformation it is Slk = hk/{2m( j £g) 1 / 2 — kXz/2. 
Once again we find the correspondence between A c and 
Xz- It follows from the above reasoning that electrons in 
a given energy band described by the one-band effective 
mass approximation are not point particles, but rather 
we should attribute to them the size Xz- 

Clearly, the magnitude of Xz is essential. There is 
Xz = A c (c/ti)(mo/m5) w 0.89(mo/mg) A since, as men- 
tioned above, u w 1.3 x 10 8 cm/sec for various materi- 
als. We estimate: for GaAs (mj w 0.067to ) A z w 13 
A, for InAs (m* w 0.024m ) X z « 37 A, for InSb 
(7Bq f=s 0.014mo) Xz ~ 64 A. Thus, in contrast to vacuum 
(A c = 3.86 x 10~ 3 A) the length of ZB in semiconductors 
can be quite large and one can confine electrons to the 
dimensions of Xz using quantum wells, wires, dots, or 
external fields. For an electron in a magnetic field the 
magnetic radius is L = (h/eB) 1 / 2 w 80 A at B = 10 T, 
and it is fairly easy to achieve L < Xz for NGS. 

It may appear unreasonable to ascribe to an electron 
the size of, say, 70 A. However, it is well known that the 
values of 1 /mj and g£ in very narrow gap semiconductors 
can be 100 times larger than the corresponding values for 
free electrons. The special property of Xz = h/m^u is 
that it cumulates the largeness of 1/toq and of \ju for 
semiconductors, as compared to 1/mo and 1/c for vac- 
uum. We emphasize again that the form of the Darwin 
term in Eq.(6) corresponds to the electron size of Xz- 

The Zittcrbcwegung length Xz appears in the descrip- 
tion of interband tunneling. The tunneling in vacuum 
from the Dirac sea of electrons with negative energies 
to the empty states of electrons with positive energies 
in the presence of a linear potential V(z) = eFz has 
been treated by RQM [19,20]. The tunneling probability 
is proportional to exp(-W), where W = nm^/ehF — 
(n/2)(zo/X c ), in which z = 2m c 2 /eF is the tunneling 
length. This probability becomes appreciable at critical 
electric fields F cr such, that zo/X c w 1, which gives F cr — 
2toqC 3 /eh ps 2.6 xlO 16 V/cm. Fields of this intensity are 
not available in terrestrial conditions. For semiconduc- 
tors the corresponding probability of tunneling between 
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FIG. 1: Energy- wave vector dependence in the forbidden gap 
of InAs. Various symbols show experimental data of Parker 
and Mead [22] for five InAs samples, the solid line is theoret- 
ical fit using Eq.(ll). The determined parameters are Xz = 
41.5 A and u = 1.33x10 s cm/sec. 

the light-hole and the conduction bands is described 

3/2 ^1/2 

by the exponential factor W — ir£ g m /2\/2ehF = 
(n/2)(zo/Xz), where z — £ g /eF is the corresponding 
tunneling length [21]. Now the critical field, at which 

zq/Xz »1, is Fz = £^ 2 m*^ 2 /y/2eh. This is of the order 
of 10 5 V/cm, which agrees with the field intensities used 
in semiconductor tunnel diodes. 

In fact, Xz can be measured directly. We write Eq.(l) 
in the form 

e = ±hu (x z 2 + k 2 ) 1/2 . (11) 

For k 2 > this formula describes the conduction and 
the light-hole bands. But for imaginary values of k there 
is k 2 < and Eq.(ll) describes the dispersion in the 
energy gap. This region is classically forbidden but it 
can become accessible through quantum tunneling. Fig- 
ure 1 shows the data for the dispersion in the gap of 
InAs, obtained by Parker and Mead [22] from tunneling 
experiments with double Schottky barriers. The solid 
line indicates the fit using Eq.(ll). The value of Xz is 
determined directly by fco for which the energy is zero: 
X z 2 = feg. The fit gives X z ~ 41.5 A and u w 1.33 x 10 8 
cm/sec, in good agreement with the above estimation for 
InAs. Similar data for GaAs give Xz between 10 A [23] 
and 13 A [24], again in good agreement with the above 
estimation. 

Our results should lead to two categories of observable 
effects. The first is related to the final electron size. For 
example, electron-electron interaction at short distances 
or electron interaction with short wavelength phonons 
should be affected by the final size (see [25]). The ef- 
fective Darwin term belongs to this group. The second 
category is associated with the question: what happens 
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when electrons are confined to dimensions smaller than 
Xz ? It is often stated in RQM that "a measurement of 
the position of a particle, such as an electron, if carried 
out with greater precision than the Compton wavelength, 
would lead to pair production" [26]. This phenomenon 
is related to our reasoning behind Eq.(3) that for such 
an electron the uncertainty of energy is larger than the 
gap. It is clear, however, that the pairs created in this 
way can only be virtual, otherwise their recombination 
would lead to the production of energy out of nothing. It 
appears that the virtual carriers could be observed, for 
example, in screening effects. 

In our theoretical treatment we used two simplifica- 
tions. The first is the use of the three-level model of band 
structure. This model is accurate near the band edges of 
NGS and it is valid up to the inflection point on the E(k) 
curve [27]. Secondly, we assumed that p x — p y = in 
the calculation of ZB and of the smearing. This was done 
for technical reasons, the difficulty being that the degen- 
eracy of the valence Tg bands couples light and heavy 
holes. However, since the conduction band is spherical 
[see Eq.(l)], it is clear that our results for the direction 
z are approximately valid for the other two directions as 
well. The isotropy of the Darwin term in Eq.(6) confirms 
this conclusion. Our simplifying assumption A>£ 9 , al- 
though not quite justified for some materials, preserves 
the essential spin properties of electrons in NGS. We add 
that another group of NGS, namely the lead salts PbTc, 



PbSe and PbS, is characterized by the 4x4 band Hamil- 
tonian closely resembling the Dirac Hamiltonian [2, 28]. 
The only complication is that the bands are ellipsoidal, 
which means that we deal with two maximum velocities, 
one for the x and y and another one for the z direction. 

To summarize, we used the band structure of nar- 
row gap semiconductors to show that, even in absence 
of external fields, semiconductor electrons experience the 
Zitterbewegung (trembling motion) with the character- 
istic frequency of loz = £g/h and the amplitude of 
Xz = h/niQU. If the electrons are described by two- 
component spinors for a given energy band (which is usu- 
ally the case) there is no ZB but the electrons should be 
treated as objects of size Xz- The magnitude of Xz in 
NGS can be as large as 70 A. Observable consequences 
of these predictions are divided into those related to the 
final electron size and those resulting from the electron 
confinement to dimensions smaller than Xz- Such obser- 
vations would lead to deeper understanding of electrons 
in solids, but they would also shed light on some still un- 
observed predictions of relativistic quantum mechanics 
and quantum field theory for electrons in vacuum. 
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